Combined Effect of Bond- and Potential-Disorder in Half-Doped Manganites 
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We analyze the effects of both bond- and potential-disorder in the vicinity of a first-order metal 
insulator transition in a two-band model for manganites using a real-space Monte Carlo method. 
Our results reveal a novel charge-ordered state coexisting with spin-glass behavior. We provide the 
basis for understanding the phase diagrams of half-doped manganites, and contrast the effects of 
bond- and potential-disorder and the combination of both. 



PACS numbers: 71.10.-w, 71.70.Ej, 75.47.Lx 

The perovskite manganites Ri_ x A^.Mn03 (R = rare 
earth, A' = alkaline earth) have received attention from 
the condensed matter community, largely due to their 
colossal magnetoresistance (CMR) effect JT], 0]. In these 
materials different ordering tendencies of charge, spin, 
lattice, and orbital degrees of freedom compete. Due 
to this complexity the understanding of the sensitive 
influence of quenched disorder has become one of the 
central issues in manganites' research. Based on a va- 
riety of experiments, Tomioka and Tokura have shown 
that the observed phases near half-doping can be or- 
ganized in terms of two parameters, the average ra- 
dius ta — (1 — x)rn + xr A' and the variance a 1 of 
the ionic radii Q. Since ta and a 1 control the single- 
particle bandwidth and the amount of disorder, respec- 
tively, quenched disorder is thereby identified as one key 
parameter in manganites. This is supported via a set of 
experiments on the "ordered" and "disordered" compo- 
sitions of Ro.5Ba .5Mn03, with a sequence of decreasing 
ta as R ranges from La to Y 0, H|. The groundstate 
for the " ordered" materials changes from a ferromagnetic 
metal (FM-M) to an A-type antiferromagnetic (AF) in- 
sulator near R=Pr, eventually followed by a charge and 
orbitally ordered insulator with ferromagnetically aligned 
spins along zig-zag chains, which is commonly referred to 
as the CE phase. In the "disordered" materials the AF 
phase does not exist and the Curie temperature (Tc) is 
suppressed. A spin-glass (SG) phase is observed for R 
ranging from Sm to Dy. "Disordered" Ndo.5Bao.5Mx1.O3 
undergoes an insulator to metal transition near Tc, which 
is characteristic of CMR materials. On applying an ex- 
ternal magnetic field the spin-glass state can be driven 
towards a ferromagnetic metallic phase 0, 0] • 

Theoretical studies of manganite-specific models at 
half-doping have indeed found a variety of ordered phases 
and transitions for clean systems [1, [|[ . Among them, the 
first-order phase transitions are particularly interesting, 
since these are assumed to split into two second-order 
transitions with an intermediate region of macroscopic 
phase coexistence upon including quenched disorder [l|. 



Detailed model studies to verify this assumption have 
only recently begun (Io| . However, the modeling of disor- 
der, which originates from different microscopic sources, 
has remained a matter of choice. Since the R or A' ions 
are located away from the electronically active Mn-C>2 
planes, the primary effect of ionic-size mismatch is to 
modify the Mn-O-Mn bond angles [111, [T^]. A realistic 
model ansatz for disorder should necessarily account for 
these bond angle variations. Additionally, the random 
positions of differently charged R 3+ and Ba 2+ ions results 
in an effective on-site disorder via Coulomb interactions. 
Therefore, both these sources of disorder are important 
and unavoidably occur together in doped manganites. 

In this Letter, we provide a systematic analysis for 
the phase diagrams of R .5Ba .5MnO3 by modeling 
the Mn-O-Mn bond angle variations by angle depen- 
dent hopping parameters in a two-dimensional two- 
band double-exchange model with electron-lattice and 
electron-electron interactions. We use the combination 
of a real-space Monte Carlo method and an unrestricted 
Hartree-Fock scheme. Bond-disorder alone explains the 
suppression of the Curie temperature and the existence 
of a spin-glass phase, but its effect on the charge and or- 
bital ordering is weak. Therefore, a phase is formed with 
coexisting charge and orbital order and glassiness in the 
spin degree of freedom. We identify a thermally driven 
metal to insulator transition, as well as a magnetic field 
driven spin-glass insulator to ferromagnetic metal transi- 
tion. We find that disorder does not lead to macroscopic 
phase coexistence, but rather generates inhomogencities 
on the scale of few lattice spacings. 

Specifically, we choose a two-band model with 
quenched disorder for itinerant e g electrons coupled to 
localized S = 3/2 t^g spins and to the Jahn- Teller (JT) 
lattice distortions. The inter-orbital Hubbard repulsion 
U' between the e g electrons and the AF superexchange 
Js between neighboring t<i g core spins are also included. 

Based on the results of previous analyses [l3[ we adopt 
the double-exchange limit for the Hund's rule coupling, 
leading to the Hamiltonian: 
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In Eq. (1), a, (3 are summed over the two Mn-e g orbitals 
d x 2_ y 2 (a) and d 3z 2_ r 2 (6). The operator c iQ (cj Q ) an- 
nihilates (creates) an electron at site i in orbital a with 
its spin slaved along the direction of the t2 g spin Si. t"? 
denote the hopping matrix elements between e g orbitals 
on nearest-neighbor Mn ions via the oxygen 2p orbitals, 
and hence depend on the Mn-O-Mn bond angle cj>ij . Tak- 
ing into account the pda contributions only, the hopping 
parameters between neighboring sites i and j are given 
by 0: 



*?/ = tcos*^ 

±ab — ±ba 



Hj ); t% = (i/3)c OS (0 y -) 
= +(-)(t/V3)~ 



(2) 



Here, x and y denote the spatial directions on a 
square lattice and t = 'i/A(pda) 2 is the basic en- 
ergy unit. The factors = cos(Qi/2)cos(Qj /2) + 
sin(@i/2)sin(Qj/2) e -1 (* i_ *i) are a consequence of pro- 
jecting out fermions with spins antialigncd to the core- 
spin directions. 0^ and <&i are the polar and azimuthal 
angles determining the orientation of the ti g spin S^. 

Bond-disorder arises from a non-uniform distribution 
of the bond-angles faj ; on-site disorder enters via the lo- 
cal potentials ej, for which equally probable values ±A 
are assumed. A denotes the strength of the JT coupling 
between the distortion = (Qf,Qf) and the orbital 



pseudospin rf 



where T M are the 



Pauli matrices [l[ . The spins are treated as classical unit 
vectors, | S»| = 1, and the lattice variables are considered 
in the adiabatic limit, /i denotes the chemical potential 
and the lattice stiffness K is set to 1. 

The average ionic radius ta in Ro.sBao.sMnOa is mod- 
elled by the average angle fa. The "ordered" compounds 
are represented by a uniform Mn-O-Mn angle faj = fa, 
and for the "disordered" systems the angles faj are se- 
lected from a binary distribution with mean fa and vari- 
ance 6<p. Since the amount of disorder depends on the 
difference in ionic radii of R and Ba ions, fa and 8<j> are 
not independent. This is reflected in the similar behav- 
ior of "ordered" and "disordered" Lao.sBan sMn03, due 
to the similar ionic radii of La and Ba [||. We there- 
fore set 5(j) = (f>Q a — fa, where fa a = 175° is the average 
bond angle in Lao.sBao.sMnOa fa is assumed to 

decrease from 175° to 164°, as R changes from La to Y. 
Our parameter choice, A = 1.4, U' = 6 and Js — 0.08 is 
guided by the observed scales for the Neel temperature 
in CaMn03 and the transport gap in LaMnOs [H, [l7| • 

Certain limits of the model Hamiltonian Eq. (1) have 
been analyzed already. The effects of large U' have 
been addressed within mean-field theories for a restricted 
choice of magnetic phases 0, [l9[ . In the absence of the 
Hubbard term, various ordered phases were discovered at 
half-doping using an exact diagonalization (ED) based 
Monte-Carlo method [|[. The effect of on-site disorder 
near a first-order phase boundary was reported recently 
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FIG. 1: (Color online) Temperature dependence of (a) the 
magnetization, (b) the spin structure factor at selected mo- 
menta, and (c) the staggered charge structure factor for clean 
and bond-disordered systems, (d) Field-cooled (FC) and zero- 
FC magnetization for the bond-disordered system. 



[nil ] . But a realistic and simultaneous treatment of bond- 
and potential-disorder and the explicit treatment of an 
inter-orbital Hubbard term, while exactly retaining the 
spatial correlations has so far been lacking. 

Here, we use a combination of the unrestricted Hartree- 
Fock (HF) scheme and the travelling cluster approxima- 
tion (TCA). The TCA allows for sampling of classical 
configurations for spin and lattice variables according 
to the Boltzmann weight, and involves iterative ED of 
small clusters [2(|. The Hartree-Fock decomposition of 
the Hubbard term leads to three independent HF param- 
eters per site (rii a ), (riib) and (t^~) — {c\ b c ia ), which enter 
as additional parameters in TCA. These HF parameters 
are self-consistcntly evaluated with the annealing process 
of the TCA. Most results are on lattices with = 16 2 
sites using a travelling cluster with N c = 4 2 . 

The indicators for spin- and charge ordering are plot- 
ted in Fig. 1. Panel (a) shows the temperature depen- 
dent magnetization, defined via m 2 = ((A -1 Si) 2 ) a «, 
for clean and bond-disordered systems for fa = 170°. 
Here and below (...) a „ denotes the average over thermal 
equilibrium configurations, and additionally over realiza- 
tions of quenched disorder. Bond disorder with S(f> = 5° 
considerably reduces the Curie temperature Tq, which 
is estimated from the inflection-point in m(T). Recall 
that the strength of bond-disorder 8(f) is tied to the av- 
erage bond angle fa via 5cj) = fa^ a — fa. In the clean 
system m(T) and in particular Tq, arc not affected 
much upon varying fa. The fingerprint of CE spin or- 
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dcr is the simultaneous presence of peaks at wavevectors 
q= (0, 7r), (tt, 0) and (tt/2, ir/2), in the spin structure fac- 
tor D 5 (q)= A^ 2 Ey(SrSj)a« e-^-C"-'*). In Fig. 1(b), 
a simultaneous rise in the three components of Dg(q) is 
observed for (50 = 0°, whereas in the disordered case no 
sizable peaks appear at any q. 

For the clean system, the charge structure factor 
£> n (q) = iV- 2 Eii(^ ' »i)a» e- i€ »-< r «-'^ at q - q = 
(jt, 7r) rises sharply upon cooling (see Fig. 1(c)), indicat- 
ing the onset of staggered charge order. The temperature 
scales Tqo and Tqe are inferred from inflection points 
in the T-dependence of the relevant components of the 
charge and spin structure factors. Charge ordering is 
accompanied by the ordering of lattice and orbital vari- 
ables. The U' term enhances Tco-, which is otherwise of 
the same order as Tqe for the values of A used here @. 
Surprisingly, a clean signal for the onset of charge order- 
ing near T ~ 0.06 is found even in the bond-disordered 
system. This highlights the crucial qualitative difference 
between the effects of bond and on-site disorder. The 
latter is known to strongly suppress charge ordering ten- 
dencies [2l|. Upon comparing Fig. 1(b) and Fig. 1(c) 
it becomes clear that the bond-disorder leads to a state 
which is ordered in the orbital and the charge sector but 
disordered in the spin sector. The U' term further stabi- 
lizes the charge order in this state. 

Additional information about the spin state that 
emerges in the presence of bond disorder is obtained from 
the field-cooled (FC) and zero-field-cooled (ZFC) mag- 
netizations, which we calculate by including a Zeeman 
term —h^i^ii m the Hamiltonian. While the FC and 
ZFC magnetizations are indistinguishable for clean sys- 
tem, their difference at low T may serve as an indicator 
for a spin-glass character for bond-disordered system (see 
Fig. 1(d)) [22|. The temperature for which the FC and 
ZFC results begin to differ, provides an estimate for the 
spin-glass crossover temperature T g . Such a state with 
charge and orbital order but glassiness in the spin sector 
was recently reported for single layered manganites [23[ . 

Fig. 2(a) shows the phase diagram in the T— O plane, 
for the clean system. For 0o > 170°, the system under- 
goes a paramagnet (PM) to fcrromagnet transition upon 
cooling. For 0o < 169°, charge and orbital degrees of 
freedom order at low temperatures, followed by a transi- 
tion from a PM to a CE spin state near T ~ 0.02. The 
0o-driven transition at low T results from the reduction 
in bandwidth upon decreasing 0o- Thereby the effective 
A and Js are enhanced leading to a first-order transition 
towards a charge and orbital ordered (CO-OO) CE state 
by opening a gap in the spectrum, at a critical value 0q . 

increases from 155° for U' = 0, to 169° for U' = 6. 
While charge and orbital order is only weakly affected 
by bond disorder, the spin degree of freedom reacts more 
sensitively (see Fig. 2(b)). For O < 169°, glassiness 
in the spin sector is induced (see Fig. 1(d)), and for 
172° > 0o > 170°, Tc is strongly suppressed. 

Some of the key features of the experimental phase di- 
agram for the "disordered" materials are not captured 
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FIG. 2: (Color online) T—<f>o phase diagrams at A = 1.4, Js = 
0.08, and U' = 6, (a) 50 = 0°, A = 0, (b) 50 = 175° - <f> , 
A = 0, and (c) <50 = 175° — 0o, A = 0.4. Symbols are data 
points and dotted lines are guides to eye. FM (PM) denotes 
a ferromagnetic (paramagnetic) state; M (I) indicates metal- 
lic (insulating) character; CO (OO) refers to charge (orbital) 
order and SG marks the spin-glass state, (d) T-dependent re- 
sistivity p for increasing strength of on-site disorder, in units 
of h/ixe 2 . p(T) for varying magnetic field for A = 0.4, and (e) 
0o = 172°, (f) 0o = 168°. 



by bond disorder alone, i.e., there is no charge disor- 
dered SG state and a FM-M to PM-I transition exists 
only in a very narrow 0o window. Therefore, we include 
additional on-site disorder. Charge order is lost for finite 
on-site disorder, while the spin-glass phase persists (see 
Fig. 2(c)). The ground state becomes an unsaturated 
FM for 170° > 0o > 173° and T c is further reduced. 
The phase diagram in Fig. 2(c) compares very well with 
the experimental phase diagram of half-doped mangan- 
ites @. 

Fig. 2(d) shows the effect of on-site disorder on the 
resistivity p of a bond-disordered system with a ferro- 
magnetic metallic ground state, p is approximated by 
the inverse of cr(u> m i n ), where lOmin = Wt/N ~ 0.04t is 
the lowest reliable energy scale for calculations of the op- 
tical conductivity <j(uj) on our 16 2 system (24|. Insulating 
or metallic character is determined from the sign of the 
slope of p(T). Increasing the strength A of the on-site 
disorder leads to a reduction in Tc and an increase in the 
resistivity, and an FM-M to PM-I transition is observed 
for A = 0.3. 

If the ground state is ferromagnetic, the resistivity in 
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FIG. 3: Monte Carlo snap-shots for (f> = 168° on a 16 x 16 
lattice. Column-(a): 5<f) — 0°, A = 0, column-(b): 5<j> — 7°, 
A = and column-(c): 8cj> = 0°, A = 0.4. The top row 
shows the charge density m, grayscale covering the range from 
0.2(white) to 0.8 (black). The bottom row plots the spin 
component S", grayscale from -1 (white) to 1 (black). Results 
are shown for specific realizations of disorder at T = 0.005. 

an external field drops near Tc (see Fig. 2(e)), simi- 
lar to the experiments in Ndo.sBao.sMnOa. A magnetic 
field reduces the resistivity in the SG phase at low T; 
additionally dp/dT changes sign, indicating an insulator 
to metal transition (see Fig. 2(f)). Such transitions are 
indeed observed in (Smo.3Gdo.7)o.55Sro.45Mn03 [H, 0|- 

Fig. 3 shows real-space patterns for the charge rii and 
the spin Sf variables for clean (column a), bond disor- 
dered (column b), and on-site disordered (column c) sys- 
tems. The groundstate in the clean case is a CE phase 
with charge and orbital order. A checkerboard pattern 
for the charge density and a zig-zag FM chain struc- 
ture for spin variables characterize this phase Q. Bond 



disorder does not affect the charge ordering, while the 
spins S become disordered leading to a nontrivial state 
with simultaneous charge order and spin glassincss. On- 
site disorder spoils the charge ordering and the magnetic 
order. The spatial patterns reveal inhomogeneities on 
the scale of a few lattice spacings. Since charge order 
can exist in the PM phase at elevated temperatures (see 
Fig. 2(a), (b)), its coexistence with spin-glass behavior 
is not surprising. The different effects of bond disorder 
on charge and spin degrees of freedom is, however, non- 
trivial. On the one hand, a decreasing bond-angle leads 
to a lowering of the fermionic kinetic energy across that 
bond and thereby reduces the tendency towards double 
exchange ferromagnetism. But more importantly, the 
zig-zag structure of the spin-aligned chains in the CE 
phase is disrupted by random weak bonds. The charge 
variables instead vary only weakly with bond disorder, 
which explains the stability of the charge ordered state. 

In summary, a two-band double-exchange model with 
electron-lattice and electron-electron interactions with 
bond and on-site disorder provides the basis for an over- 
all understanding of the experimental phase diagram of 
Ro.5Bao.5Mn03 Q. Our results reveal the existence of 
a non-trivial phase, which is ordered in the charge vari- 
ables but glassy in the spin variables. The bond disorder 
is important for the spin-glass behavior, and likely to 
explain the coexistence of charge order and spin glassi- 
ness in single- layered manganitcs 23] . On-site disorder is 
crucial for a description of the thermally driven metal to 
insulator transitions. Each of the two types of disorder 
explain selected features in manganites, but only a com- 
bination of both can describe most of the experimentally 
observed phenomena. 

We acknowledge support by the Deutsche Forschungs- 
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owulf Cluster at HRI, Allahabad (India). 
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